arXiv:1505.02770vl [astro-ph.CO] 11 May 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 


Printed 13 May 2015 (MN MbX style file v2.2) 


Semi-Analytic Galaxy Formation in coupled Dark Energy 
Cosmologies 

Fabio Fontanot^*, Marco Baldi^’^’^, Volker Springel^’® and Davide Bianchi^ 

^ INAF - Astronomical Obserx’atory of Trieste, via G.B. Tiepolo 11,1-34143 Trieste, Italy 

^ Dipartimento di Fisica e Astronomia, Alma Mater Stitdiorum Universit‘a di Bologna, viale Berti Pichat, 6/2, 1-40127 Bologna, Italy 
^ INAF - Osservatorio Astronomico di Bologna, via Ranzani 1,1-40127 Bologna, Italy 
INFN - Sezione di Bologna, viale Berti Pichat 6/2,1-40127 Bologna, Italy 
® HITS-Heidelberger Institutfiir Theoretische Studien, Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany 
® Zentrum fiir Astronomie der Universitdt Heidelberg, ARl, Mbnchhofstrasse 12-14, 69120 Heidelberg, Germany 
^ Institute of Cosmology ami Gravitation, University of Portsmouth, Portsmouth POl 3FX 


Accepted ... Received ... 


ABSTRACT 

Among the possible alternatives to the standard cosmological model (ACDM ), coupled Dark 
Energy models postulate that Dark Energy (DE), seen as a dynamical scalar field, may in¬ 
teract with Dark Matter (DM), giving rise to a “fifth-force”, felt by DM particles only. In 
this paper, we study the impact of these cosmologies on the statistical properties of galaxy 
populations by combining high-resolution numerical simulations with semi-analytic models 
(SAM) of galaxy formation and evolution. New features have been implemented in the ref¬ 
erence SAM in order to have it run self-consistently and calibrated on these cosmological 
simulations. They include an appropriate modification of the mass temperature relation and 
of the baryon fraction in DM haloes, due to the different virial scalings and to the gravitational 
bias, respectively. Our results show that the predictions of our coupled-DE SAM do not differ 
significantly from theoretical predictions obtained with standard SAMs applied to a reference 
ACDM simulation, implying that the statistical properties of galaxies provide only a weak 
probe for these alternative cosmological models. On the other hand, we show that both galaxy 
bias and the galaxy pairwise velocity distribution are sensitive to coupled DE models: this 
implies that these probes might be successfully applied to disentangle among quintessence, 
/(i?)-Gravity and coupled DE models. 
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1 INTRODUCTION 

Dark Energy (DE) represents a critical unknown for the concor¬ 
dance cosmological model of our Universe, which emerged as the 
result of a decade long effor t in the determination of the key co smo- 
logical parameters (see e.g IPlanck Collaboration X^l2014h . The 
easiest description for this mysterious contributor to the present 
energy density, which accounts for ~ 70 per cent, is a classical 
cosmological constant A, i.e. an homogeneous and static energy 
density filling the whole Universe at all epochs. This simple model 
(ACDM hereafter), while indeed able to explain the vast majority 
of the observed propertie s of the Univers e, bears a number of theo¬ 
retical problems (see e. g. IWeinberj 1 98^ for a review), mostly due 
to the level of “fine-tuning” required to accommodate for the small 
value of A at the present epoch. A number of alternative DE mod¬ 
els have thus been proposed in the literature trying to explain the 
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origin of the accelerated expansion; these range from scalar field 
theories (i.e. quintesse nce), to modifica tions of the equation of gen¬ 
eral relativity (see e.g. lAmendolal201^ and references therein for a 
comprehensive review). In order to provide the observational con¬ 
straints needed to disentangle between those different scenarios, 
wide galaxy survey s are currently unde r advanced planning, like 
the Euclid mission dUaureiis et alJllOllh . which relies on a combi¬ 
nation of weak lensing measurements (based on precision imaging) 
and clustering analysis (from slitless spectroscopy), which will al¬ 
low to study, at the same time, both the evolution of the equation of 
state of DE and the growth function. 

Since galaxies provide the privileged tracers of cosmic evolu¬ 
tion, it is crucial, for the success of these missions, to understand 
the interplay between the physical processes responsible for galaxy 
formation and evolution and the assembly of the large-scale struc¬ 
ture. The former issue has been explored, for a number of alterna¬ 
tive cosmologies, with the help of numerical methods that follow 
the non-linear evolution of virialized structures (i.e. high resolu¬ 
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tion A ^ -body simulations, see e.g. iGrossi & Springeill2009l : iBaldil 
l2012bl : IPuchwein et alj|?013l and references herein). In the bary- 
onic sector, semi-analytic models (SAMs) employ simplified an¬ 
alytic prescriptions to model the relevant processes acting on the 
baryonic gas (and their interplay), and to study the evolution of 
galaxy components as a function of their physical properties, red- 
shift and environment. This approach has been shown to correctly 
reproduce a number of observational data, a nd the predictions of 
differ ent models are consistent in many cases l lpontanot et al.l200^ . 
2012h . even if some tension with the data still remains (see e.g. 


McCarthy et al.l200^:lBovlan-Kolchin et alj20lalWeinmann et al.1 


2012h . The analytic prescriptions embedded into SAMs are phys¬ 
ically grounded and observationally motivated but involve numer¬ 
ous parameters, which are constrained by requiring the model to 
reproduce a well defined sample of (typically) low-z observations. 
Thi s approach thus harbou rs a significant level of degeneracies (see 
e.g lHenrigues et aDl2009l) . which are also related to the fact that 
different authors adopt different approximations for key physical 
mechanisms. In the context of future space missions, it is therefore 
of fundamental importance to characterize the impact of alternative 
DE models on the predicted properties of galaxy populations, in or¬ 
der to devise observational tests that can safely distinguish different 
cosmological models. 

Our group was the first to study the implications of non¬ 
standard cosmological scenarios on the relevant properties of 
galaxies in a large cosmological volume, and to fully quantify their 
impact on the statistical properties of galaxy populations as pre¬ 
dicted by SAMs. We focus on the amplitude of the expected modi¬ 
fications in the galaxy stellar mass function, on the cosmic star for¬ 
mation rate and on the 2-point correlation function, but we also con¬ 
sider higher order statistics like galaxy bias and the galaxy pairwise 
velocity distribution. For the cosmologies we tested so far, we con¬ 
cluded that galaxy properties alone are usually inefficient to con¬ 
strain cosmological models beyond ACDM , but whenever they are 
combined with suitable information on the underlying Dark Matter 
(DM) distribution, it is possible to devise statistical tests able to 
disentangle these alternative cosmological scenarios from the stan¬ 
dard model and among themselves. This is the four th paper on our 
series : we already considered Early D ark Energy dFon tanqt et ^ 
I 2 OI 2 L hereafter Paper I), /(i?)-Gravity jpontanot et alj201^. here¬ 
after Paper II) and massive neutrino cosmologies 1 Fontanot et al.l 
I 2 OI 5 L hereafter Paper III). In this work, we consider a new class of 
cosmologies, namely the so-called coupled DE m odels (cDE here¬ 
after, see e.g. IWettericlillI 995 I : lAmendolallTOOOh . which are based 
on the dynamical evolution of a classical scalar field 4> that plays 
the role of DE, and interacts directly with the Cold Dark Matter 
(CDM) fluid. As a consequence of the exchange of energy and 
momentum during cosmic evolution, the DE scalar field thus me¬ 
diates a “fifth-force” between Dark Matter particles, leading to a 
modification of the gravitational growth process (both at linear and 
non-linear scales) and to a different evolution of cosmic large-scale 
structure. N-body simulations of cDE cosmologies h ave been per¬ 
formed by various groups in the last decade (see e.g. Maccib et al.l 


l2004lBaldi et al.l2()ld:Ei & Barrowl20lII:ICarlesi et alj2014h and 

represent now a robust tool to investigate DE interactions in the 
non-linear regime. The particular imprints of th is class of mod¬ 
els have been studied in detail by iBaldil l l2012bh using a suite of 
W-body and hydrodynamical simulations (the CODECS projecQ), 


^ The simulations data of the CODECS project are publicly available at 
www.marcobaldi.it/CoDECS 


Table 1. Cosmological parameters for our simulations. The columns con¬ 
tain from left to right: the normalization A and the exponent a of the poten¬ 
tial, the strength of the coupling /3, as and the equation of state parameter 
ru at 2 = 0. 



A 

OL 

/3 

o' 

II 

bi 

GO 

b 

II 

0 

ACDM 

— 

— 

— 

0.809 

-1.0 

EXP003 

0.0218 

0.08 

0.15 

0.967 

-0.992 

SUGRA003 

0.0202 

2.15 

-0.15 

0.806 

-0.901 


focusing on the degeneracy of DE coupling with other cosmologi¬ 
cal parameters (like the normalization of the matter power spectrum 
as) and on the gravitational bias, i.e. the offset between the density 
fluctuations in CDM and in the baryonic components (leading to 
different baryon fraction at cluster scales). 

This paper is organized as follows. In Section[2 we introduce 
the cosmological numerical simulations and semi-analytic models 
we use in our analysis. We present the predicted galaxy properties 
and compare them among different cosmologies in Section Fi¬ 
nally, we discuss our conclusions in Section|4] 


2 MODELS 

2.1 Coupled Dark Energy cosmologies 

In this work, we consider a set of flat cosmological models includ¬ 
ing CDM, baryons, radiation and a DE scalar field 0. Among the 
range of models included in the CODECS project we focus on two 
different choices for t he scalar field po tential. We first consi der an 
expone ntial potential ( IWetterichlll988L the EXP003 model in iBaldil 

l2012bh 

!/(<(.) = ( 1 ) 


which is characterized by stable scaling solutions for the scalar 
field independently from initial conditions. In particular, in cDE 
scenarios, such a potential provides a transient early DE solutio n 
and a late time accelerated attractor (s ee e.g. Amendolal 20041) . 
We also consider a SUGRA pote ntial jBrax & Martini 1 19991 the 
SUGRA003 model in lBaldill2012bl) 


- /2 


= Acj,-°‘e 


( 2 ) 


which is typical of supersymmetric theories of gravity and implies 
a “bounce” of the DE e quation of sta te w at the cosmological “bar¬ 
rier” u) = — 1 (see e.g. lBaldill2012all . This feature has relevant im¬ 
plications for the expected number density evolution of DM haloes, 
as well as for the evo lution of the cosmological Hubble function 
and growth factor (see lBaldll2012bL for more details). 

The evolution of the main cosmological components is de¬ 
scribed by a set of dynamical equations including the interaction 
between the scalar field and CDM particles (i.e. the right-hand side 
of Equations and |4j : 


<l> + 3H4> + 


dcj) 


IGttG 


/3(0)pc 


(3) 


Pc + 3T/pc = - yi^^(</,)pc0 (4) 

where pc represents the density of CDM particles and an overdot 
denotes the (cosmic) time derivative. In the cDE models consid¬ 
ered in this work, the coupling function /3((/)), which controls the 
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interaction strengti0, is assumed to be constant. As baryons and 
radiations are always uncoupled from cj) their evolution follows the 
usual relations ph oc a~^ and pr oc a~‘^, respectively. 

At the level of linear perturbations, the interaction determines 
a modification of the growth rate due to the presence of a fifth-force 
with a strength 4/3^/3 times the standard gravitational accelera¬ 
tion acting between CDM particles, and to an additional velocity- 
dependent acceleration proportional to 0/3. These effective modifi¬ 
cations of the standard gravitational evolution affect also the non¬ 
linear dynamics of collapsed structures. For a more detailed discus¬ 
sion on the linear and non-linear properties of cDE cosmo logi es, we 
refer t he reader to lAmendoi3 l l2004l) . lBaldi et alj j2010l) and iBaldil 

( l201lh . 


2.2 Numerical simulations 


In this work, we tak e advantage o f the results of the CODECS nu¬ 
merical simulat ions llBaldill2012hll. using a modified version of the 
GADGET code dSpringel et al.l 2005h . designed to inclu de the spe¬ 
cific physical processes arising in the cDE scenario jBaldi et al.l 
I 2 OIOI 1 . In particular, we analyse the outcome of the H-CoDECS 
set, i.e. adiabatic hydrodynamical simulations on periodic boxes 
80h“^Mpc on a side, using 2 x 512® particles (correspond¬ 
ing to a mass resolution of 2.39 x 10® for CDM and 


4.78 X 10^ h“®MQ for baryons). An entro py-conserving formu¬ 
lation of smoothed particle hydrodynamics dSoringel & Hemguisll 
I 2 OO 2 L SPH) has been used to estimate hydrodynamical forces act¬ 
ing on gas particles, and no additional radiative processes (gas 
cooling, star formation, feedbacks) have been included. Initial con¬ 
ditions for all runs were generated using N-GENIC by displac¬ 
ing particles from a homogeneous glass distribution imposing the 
same amplitude of the initial power spectrum at the last scatter¬ 
ing surface, the same phases and mode amplitudes, to ensure a 
similar realisation of the large scale structure and to allow an 
object-hy-object comparison. For all simulations, a flat cosmolog¬ 
ical model has been assumed with z = 0 cosmological parame¬ 
ters consistent wit h the 7th year results o f the Wilkinson Microwave 
Anisotropy Probe jKomatsu et al.ll2oTll . WMAP7), i.e. density pa¬ 
rameters flcDM = 0.226, Ddb = 0.729 and flbar = 0.0451 (for 
CDM, DE and baryons respectively), Hubble parameter h = 0.703 
and Gaussian density fluctuations with a scale-invariant primordial 
power spectrum with spectral index n = 0.966. The common nor¬ 
malisation of perturbations at CMB and the different growth factors 
for cDE runs imply that the EXP003 run has a different amplitude 
of density perturbation at every redshift and a different cts at 2 = 0, 
as listed in Table[T] On the other hand, the SUGRA003 run has, by 
construction, the same linear normalisation as ACDM both at CMB 
and at present, with the final results that its crs value is very similar 
to the ACDM run. 

For each run, 69 snapshots were storec0; the corresponding 
group catalogues were generated using a Friend-of-Friend algo¬ 
rithm with a linking length of 0.2 (in mean particle separation 
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Figure 1. Redshift evolution of the normalization of the virial scaling re¬ 
lations for DM haloes (k, upper panel) and baryon fractions (/bar. lower 
panel). In each panel the black crosses, red diamonds and green stars refer 
to the ACDM , EXP003 and SUGRA003 runs respectively. Solid, dashed 
and dot-dashed lines represent the 3-rd order polynomial best fits for each 
cosmology as indicated in the legends. 


units), and gravitati onally bound substru ctures have been identi¬ 
fied using SUBFIND dSoringel et alJl200lh (only subhaloes that re¬ 
tain at least 20 particles after the gravitational unbinding procedure 
were considered). We then u sed the subhalo catal ogues to define 
the merger tree histories as in lSpringel et alj ( l2005h . 


2.3 Semi-Analytic Models 


In this paper, we consider the same SAM suite we used in the pre¬ 
vious papers of the series. This includes three different versions 
of the L-GALAXIES semi-analytic model, n a mely those presented 


Croton et alj (l2006h . IPe Lucia & Blaiz^ (l2007h and IGuo et akl 
(I 2 OI ll) . All these models were run by construction on Millennium¬ 
like merger tree as defined in the previous section and they repre¬ 
sent a coherent set of modelfl suitable to study the intra-model 
variance due to different assumptions made in the modelling of 
the key physical processes. The free parameters usually associated 
with this approach have been calibrated (for a ACDM cosmologi¬ 
cal model), by comparing model predictions to a well defined set 
of observational constraints (mainly at low-redshift). In order to 
highlight the differences in galaxy properties due to the different 


^ In particular, the .sign of 0/3(0) is related to the energy-momentum flow 
between the two components, such as negative (positive) values of this 
quantity conespond to a transfer from DE to CDM (from CDM to DE) 
and to an increase (decrease) of the DM particle mass. 

® At variance with our previous work, the snapshot list is not the same as 
in the Millennium simulation; the CODECS redshift sampling was chosen 
mostly to allow the constructio n of full light-cone s for weak lensing and 
CMB lensing purposes (see e.g. lGiocoli et alJlOlSl and references herein) 


Thus avoiding any additiona l noise in the predic tions due to different 
definitions of DM merger trees iKnebe & et al.l20T^ . 

® These three mod els mark the h i storica l evolu tion of the c ode o rig¬ 
inally developed by ISpri ngel et alj l2005ll : from ICroton et alj l200d) to 
IPe Lucia & Blaizog j2007h the main changes involve the treatment of dy¬ 
namical friction and merger times, the initial m ass function (from Salpeter 
to Chabrier) and the dust modelling; while from lPe Lucia & Blaizoll fcOOTh 
to lGuo et alj |201 ih the modelling of supemovae feedback, the treatment of 
satellite galaxy evolution, tidal stripping and mergers were added. 
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Figure 2. Predicted galaxy properties for different coupled Dark 


(light grey points refer to the compilation from iFontanoteljf 


model and the corresponding mass function from thelGuTeT^ 


Energy cosmological scenarios. Upper panels: redshift evolution of the stellar mass function 


200S)). The lower row shows the ratio between the mass function in a given cosmological 


201 ih model in the ACDM run. Lower panel: Cosmic star formation rate density (light grey 


points refer to the compilation from lHopkinj i20o3l l. In each panel the solid black, dashed red and dot-dashed green li nes refer to SAM predictions in ACDM 
, EXP 003 and SUGRA003 cos mologies respectively. Dark grey areas mark the distribution of the predictions from the lGuo et al.l l201lh . lDe Lucia & BlaizoJ 
i2007h and ICroton et alj 120061) SAMs applied to the same ACDM run. 


underlying cosmology, we choose not to recalibrate the models, 
thus holding the role of other astrophysical processes fixed. This 
implies that the models are optimally calibrated to reproduce ob¬ 
servations only for the ACDM run, where we expect the scatter in 
their predictions to be representative. 


a modified version of the I Guo et al.l ll201lh modeQrun on the cDE 
boxes. The main changes with respect to the ACDM SAM ver¬ 
sion include the following new features. First, the code handles a 


In the next section, we present the predictions obtained from 


° For the s ake of simplicity, in the following, we still refer to our modified 
code as the lGuo et alj 1201 ih model 
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user-generated Hubble function (tabulated in an external file) ex¬ 
tracted from the corresponding numerical realization of the cDE 
cosmology under analysis. In Paper II, we showed that, for cos¬ 
mologies introducing a fifth force, the viri al scaling relations devi- 
ate from the ACDM expectations (see also lArnold et al.ir2014h . i.e. 
the relation between total DM mass inside a sphere with interior 
mean density 200 times the critical density at a given redshift and 
the one-dimensional velocity dispersion inside the same radiufl is 
modified. As in Paper II, we find that the actual relations in cDE 
cosmologies (defined using > 10^^ Mq DM haloes) are offset ver¬ 
sions of ACDM ones (see e.g. Fig. 1 in Paper II). Nonetheless, at 
variance with the constant shift for unscreened haloes we found in 
the /(/?)-Gravity case, the offset for cDE cosmologies is redshift 
dependent, and smaller in amplitude (roughly corresponding to a 1 
percent variation at most). We show the different values for the nor¬ 
malization K of the fitted virial relations in Fig.[T](upper panel). As 
in Paper II, we then apply a correction to the virial scalings assumed 
in L-GALAXIES , corresponding to the ratio between the virial scal¬ 
ings in the desired cosmology and the corresponding virial scalings 
in the ACDM run with the same initial conditions. We then model 
the evolution of the offsets as a 3-rd order polynomial as a function 
of redshift. 


A specific feature of cDE co smologies is the different baryon 
fraction (/bar) l lBaldi et alj2010h in massive haloes, due to the dif¬ 
ferent forces felt by the baryons and DM particles. We estimate 
/bar for > DM haloes in our simulations and show the 

redshift evolutio n of this quantity i n Fig. [T] (bottom panel). Our re¬ 
sults confirm the lBaldi et alj j2010h results; we also stress that /bar 
evolves as a function of redshift, and it does not show any DM 
mass dependence at fixed cosmic epoch. Typical differences with 
respect to the ACDM run are smaller than 3 percent at z < 2. 
In L-GALAXIES , /bar mainly regulates the infall of pristine gas 
when DM haloes grow by accretion from the surrounding field and 
it is usually treated as a redshift independent free parameter. In our 
modified code, we still keep the baryon fraction as a free parameter, 
but we require it to scale, at a given redshift, as the ratio between 
the corresponding /bar in the cDE and ACDM boxes. Also in this 
case, we model the baryon fraction evolution in different cosmolo¬ 
gies as a 3-rd order polynomial in redshift. It is worth stressing, 
that the /bar shown in Fig. [T] refers to the ratio between the mass 
in baryons and the DM mass, yielding a definition closer to the 
quantities actually used in L-GALAXIES . We also consider the al¬ 
ternative definition of /bar as the ratio between the baryonic mass 
and the total mass in the halo: for all cosmological models, the dif¬ 
ferent definition mainly changes the normalization of /bar ( 2 ), but 
not its evolution. Therefore, as our modifications to L-GALAXIES 
involve only fractional quantities, our results are insensitive to the 
/bar definition employed. 


It is also worth stressing that these two features we included 
in the SAM have a limited (but not negligible) impact on model 
predictions, due to the small offsets with respect to the ACDM re¬ 
alization and that most of the differences we will discuss in the 
following are triggered by differences in merger trees and cosmic 
growth history of the large scale structure realized in the cosmolog¬ 
ical volumes under analysis. 



-1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 

Log(r / h"‘ Mpc) 


Figure 3. Redshift evolution of galaxy bias. In each panel, only model 
galaxies with M* > IO^Mq have been considered. Line types, colours 
and shaded areas have the same meaning as in Figure|2] 


3 RESULTS & DISCUSSION 


In Figure [J] we present a selection of the statistical properties of 
M* > 10® Mq galaxies, as predicted by the different SAMs. In 
the upper and lower panels, we show the redshift evolution of the 
galaxy stellar mass function and cosmic star formation rate, respec¬ 
tively. The model predictions have been convolved with a lognor¬ 
mal error distribution with amplitude 0.25 (0.3) for stellar masses 
(star formation rates) to account for the typi cal observational erro r 
in the estimate of these physical quantities (|Fbntanot_et^|200^. 
Stellar masses and star formation rates in the Croton et al.l 1 20061) 
model have been converted from Salpeter to Chabrier IMF ap¬ 
plying a constant shift of 0.25 and 0.176 dex, respecti vely. In all 
panels , solid black lines refer to the predictions of the IGuo et al.l 
1 I 2 OI ih model for the ACDM simulation. We also consider the pre¬ 
dictions of the other two SAMs applied to the same simulation to 
estimate the scatter in the predictions of different SAMs applied 
to the same cosmological box (shaded areas): we stress again that 
the main source of this scatter lies in the different treatment of 
the key physical mechanisms driving galaxy evolution in the three 
mo dels. In all panels , we also show the predictions for the modi- 
fied lGuo etallooilh model in the EXP003 and SUGRA003 runs, 
with red dashed and green dot-dashed lines, respectively. As far as 
galaxy properties are concerned, the deviations of the model pre¬ 
dictions in these alternative cosmological scenarios with respect to 
the reference ACDM run are quite small and comparable to the 
intra-model variance at fixed ACDM cosmology. This provides an 
a-posteriori justification for our choice not to recalibrate the model 
parameters, as the properties of the overall galaxy populations are 
consistent among the different cosmological runs. There is a slight 


^ lEvrard et al.l J2008h showed that haloes in ACDM cosmology are ex¬ 
pected to follow the theoretical scaling (T 200 oc [h{z)M20o{z)]^ 1 ^ 
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V||(1+z)/H(z) [Mpc/h] 


Figure 4. Pairwise galaxy velocity distribution along the line of sight for the cDE cosmologies considered in this paper at four different redshifts. Different 
cosmological models are marked by different linetypes, shading and colours as in Figure|2] The upper left caption in each panel states the different combination 
for values of the galaxy separation (r^,r||), perpendicular and parallel to the line of sight, respectively; in all panels, the conformal Hubble function l-L = aH 
has been used to rescale velocities to comoving distances. 


tendency for the deviations to grow with redshift, due to the differ¬ 
ent CT8 value at z = 0. At variance with the /(7?)-Gravity cosmol¬ 
ogy we studied in Paper II, we find no relevant dependence of these 
deviations on stellar mass. In Paper II, we interpret this effect as be¬ 
ing due to the different virial scalings in the /(i?)-Gravity model, 
which directly affect the modeling of AGN feedback. Here the de¬ 
viations from the ACDM virial scalings are definitely smaller than 
in the /(7?)-Gravity case, and the lack of a stellar mass dependence 
in the ratio between the mass function in a given cosmology and in 


the ACDM run, clearly supports our conclusion that this effect is 
negligible for the cDE cosmologies we consider in this paper. All 
the differences seen in Fig. are driven by the different merger 
trees statistics associated with the different cosmologies. Overall, 
the deviations from the ACDM mass function remain within 0.2 
dex at most mass scales and redshifts. 

In our previous work, we considered two additional cosmolog¬ 
ical tests, namely galaxy bias and the pairwise velocity distribution, 
and we discussed their efficiency in disentangling between ACDM 
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and other cosmological scenarios. The discriminating power of 
these observables is mainly driven by the combination of galaxy 
populations statistics with information on the distribution of the 
DM in the underlying large scale structure. In Figure [3] we show 
the galaxy bias estimates for the cosmological models considered 
in this work based on the ratio between the auto-correlation func¬ 
tion of galaxies in real space ^gai and the corresponding ^cdm com¬ 
puted for a randomly selected subsample of 1 percent of the CDM 
particles in the cosmological box. Only the EXP003 model shows 
a deviation from the ACDM model larger than the SAM variance 
in the standard cosmology: in particular this holds at 2 : > 1 and 
for scales larger than ~1 /i“^Mpc. The different as in the EXP003 
run contributes to this deviation, but it cannot totally account for it 
at sca l es sma ller than ~10 /i“^Mpc and for 2 > 1, as shown by 
ISaldil (l2012bl . their Figure 5). Also the SUGRA003 model shows 
clear deviations from the corresponding ACDM run, but those are 
of the same order as the intra-SAM variance, thus limiting the effi¬ 
ciency of this estimator in breaking the degeneracies between dif¬ 
ferent cosmologies. 

Moreover, in Figure |4] we present the redshift evolution for 
the pairwise galaxy velocity distribution along the line of sight 
V{v \\, r||, rx), where we consider fixed components of galaxy sep- 
aration parallel (rn ) and perpendicular (r±) to the line of sight 
( IScoccimarrol2004l . see e.g.). This quantity is a reliable indicator of 
the assembly of the large scale structure, as it traces the anisotropy 
of redshift-space correlation functions and it is strongly sensitive 
to the abundance of massive haloes. In Figure|4]we adopt the same 
reference separations we choose in our previous work (i.e. I and 
15 /i“^Mpc), despite the smaller box size of the H-CoDECS runs 
(i.e. cubic boxes with 80/i~^Mpc sides). Also in this plot, only 
galaxies with M* > 10 ®Mq have been considered; furthermore, 
the velocities have been rescaled using the conformal Hubble func¬ 
tion H = aH in order for the distribution to represents the statis¬ 
tical displacement of galaxy pairs from real to redshift space. We 
assume the usual convention that the pairwise velocity is positive 
when galaxies are receding and negative when they are approach¬ 
ing. This analysis leads us to similar conclusions with respect to 
Figure [3l the EXP003 model clearly show a different velocity dis¬ 
tribution with respect to predictions in the ACDM run (due in part, 
but not completely, to the increase in as), while the SUGRA003 
results are virtually indistinguishable from those obtained for the 
standard cosmology. 


4 CONCLUSIONS 

In this paper, we present an updated version of the L-GALAXIES 
semi-analytic model, designed to run self-consistently on the 
Codecs sui te of numeric al simulation for coupled Dark Energy 
cosmologies ( lBaldill2012bll . The main modifications with respect 
to the ACDM version of the code include: (a) an implementation 
of an user-defined Hubble function; (b) modified DM halo virial 
scalings (which impact the mass-temperature relation); (c) modi¬ 
fied baryon fractions accounting for the gravitational bias. Item (a) 
was already introduced in Paper II; item (b) has been modified for 
this project by allowing a redshift dependence for the virial scal¬ 
ing normalization (with respect to the ACDM expectations); item 
(c) has been introduced in this work. All these new features of the 
model have been directly calibrated on the cDE simulations under 
consideration, by comparison with the corresponding ACDM runs. 

The new code represents a step forward with respect to pre¬ 
vious versions, designed to run on Early Dark Energy (Paper I) 


and /(iZ)-Gravity cosmologies (Paper II), allowing us to present 
the first comprehensive picture of the impact signature of non¬ 
standard coupled-DE cosmologies on the statistical properties of 
galaxy populations. The weak cosmological constraints coming 
from direct comparison with model galaxy properties (alone) con¬ 
firm the robustness of the SAM predictions a gainst small varia¬ 
tions in the cos mological framework (see also IWang et aljHOOSi ; 
ICuo et ai]|201 ih . On the other hand, the modification of the cos¬ 
mological framework we consider has a substantial impact on the 
growth and assembly of the large scale structure: therefore, com¬ 
bining predicted galaxy properties with information on the distri¬ 
bution of the underlying total mass distribution, it becomes possi¬ 
ble to disentangle different cosmological scenarios. In particular, 
we focus on standard tests like galaxy bias and the galaxy pair¬ 
wise velocity distribution, and we show that coupled DE models 
can be distinguished from ACDM runs using both probes, unlike 
quintessence models (in Paper I we showed that only bias is a sen¬ 
sitive probe) and /(i?)-Gravity models (in Paper II we showed that 
only the pairwise velocity distribution is a sensitive probe). How¬ 
ever, our results also show that these cosmological tests are suffi¬ 
ciently sensitive only for a subset of coupled DE cosmologies, in 
particular for the exponential potential rur0, while the run includ¬ 
ing a supersymmetric potential corresponds to a much weaker (al¬ 
though coherent) signal. Moreover, in Paper III we showed that the 
likely existence of a massive neutrino background implies devia¬ 
tions from a pure ACDM run which go on in the opposite direction 
with respect to these results, i.e. an increased galaxy bias and a 
narrower pairwise distribution. Therefore, the inclusion of such a 
background would have the net effect of reducing the cosmological 
signal com ing from either co upled DE or /(i?)-Gravity models (as 
seen i.e. in iBaldi et alj20lil . 

Overall, the results we presented in our series of papers are 
of particular relevance for the plannin g and exploitation o f future 
wide area galaxy surveys (like Euclid, Laureiis et alj|20llh . meant 
to shade light on the true nature of DE. In fact, the wealth of data 
coming from these efforts requires careful calibration and analysis 
in order to provide a proper characterization of the large scale struc¬ 
ture evolution. Of course, a deep understanding of all the systematic 
effects, either due to galaxy formation physics or the cosmological 
parameters, is of critical importance, given the exquisite precision 
required for disentangling the different scenarios. Therefore, the 
construction of mock galaxy catalogues covering the widest range 
of proposed DE theories represents a key step. Our SAM suite, de¬ 
signed to run self-consistently on a variety of these DE models, 
covers this need and provides a tool to test the relative efficiency 
of cosmological probes. In this paper, as in our previous work, 
we focus on scales suitable for galaxy studies (i.e. stellar masses 
10®Mq < Matar < IO^^Mq), but using a moderate volume: we 
thus plan to apply our SAMs to larger cosmological volumes in 
the future and use these runs to build cosmological light cones (see 
e.g. iMerson et alj|201^ and mock galaxy catalogues resembling, 
as close as possible the expected galaxy properties in the different 
DE cosmologies. 


® It is worth stressing that the particular model tested in this paper 
(EXP003) has a relatively strong value of coupling, excluded at about 3 (t 
C.L. by the most recent CMB constrains. However, such large coupling val¬ 
ues might be still viable in the presence of a substantial (but still reasonable) 
contribution of massive neutrinos to the cosmic budget. 
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